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In paper I of this series on fluid turbulence we showed that exact resummations of the perturbative 
theory of the structure functions of velocity differences result in a finite (order by order) theory. 
These findings exclude any known perturbative mechanism for anomalous scaling of the velocity 
structure functions. In this paper we continue to build the theory of turbulence and commence 
the analysis of nonperturbative effects that form the analytic basis of anomalous scaling. Starting 
from the Navier-Stokes equations (at high Reynolds number Re) we discuss the simplest examples 
of the appearance of anomalous exponents in fluid mechanics. These examples are the nonlinear 
(four-point) Green’s function and related quantities. We show that the renormalized perturbation 
theory for these functions contains “ladder" diagrams with (convergent!) logarithmic terms that 
sum up to anomalous exponents. Using a new sum rule which is derived here we calculate the 
leading anomalous exponent and show that it is critical in a sense made precise below. This result 
opens up the possibility of multiscaling of the structure functions with the outer scale of turbulence 
as the renormalization length. This possibility will be discussed in detail in the concluding paper 
III of this series. 

PACS numbers 47.27.Gs, 47.27.Jv, 05.40.+j 


1. INTRODUCTION 

In this paper we clarify, on the basis of an analytic the¬ 
ory, how anomalous scaling appears in fluid mechanics. 
The aim of the analytic theory is to reach understand¬ 
ing on the basis of the Navier-Stokes equations. We thus 
make a clear cut distinction between the analytical ap¬ 
proach and the host of ad hoc models that were employed 
to attempt a description of anomalous scaling in turbu¬ 
lence. 

The common wisdom about anomalous scaling in tur¬ 
bulence can be found in a number of recent reviews and 
books (see, e.g. In terms of the scaling exponents, 

most attention was given to those associated with the 
structure functions of velocity differences and to the cor¬ 
relation function of the fluctuations in the rate of energy 
dissipation. In terms of the Eulerian velocity field u(r, t) 
the structure functions S n ( R ) can be defined || in terms 
of the velocity differences du(r + R, r, t): 

<5u(r + R, r, t) = [u(r + R, t) — u(r, t)] , (1.1) 

S n (R)=( |(Su(r + R,r,t)|"> (1.2) 

where the symbol (...) denotes an average over time and 
over r. The structure functions S n (R) are expected to 
exhibit scaling behavior for values of R in the “inertial 
range” L 3> R o where L and ij are the outer scale 
and the Kolmogorov dissipation scale respectively: 

S n (R) ~ CR U . (1.3) 

In the Kolmogorov theory (K41) the exponents C, n equal 
n/ 3, whereas experiments and popular belief ascribed 


anomalous values to these exponents. The energy dis¬ 
sipation rate e(r,t) is the field 


e(M) = 2 


d a up(r, t) + dpu a {r , t) 


(1.4) 


where v is the kinematic viscosity. The “centered” cor¬ 
relation function of the fluctuations in this field, I\ ee (R), 
is defined as 


Kee(R) = (£(r + R, t)e( r, t)) , (1.5) 

where e(r, t) = e(r, t) — e, and e is the mean of the dis¬ 
sipation field. It was found in experiments that K ee (R) 
decays very slowly in the inertial range, 

Kee(R) ~ R (1-6) 

with /i having a numerical value [Q in the range 0.25 ± 
0.05. It was claimed (U that the K41 theory required p, to 
vanish. Accordingly, there have been many attempts to 
construct models of turbulence (see, i.e. 
and references therein) to take ( |l.G|) into account and 
to explain how measured deviations (£„ — n/ 3) in the 
exponents of the structure functions were related to /r. 

In this paper we show that the renormalized pertur¬ 
bation theory for correlation functions that include ve¬ 
locity derivatives (to second or higher power) exhibit in 
their perturbation expansion a logarithmic dependence 
on the viscous scale 77 |jl5|,^6|. In this way the inner scale 
of turbulence appears explicitly in the analytic theory. 
The perturbative series can be resummed to obtain inte- 
grodifferential equations for some many-point objects of 
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the theory. These equations have also non-perturbative 
scale-invariant solutions that can be represented as power 
laws of rj to some exponents A. For example it will be 
shown in paper III that the correlation function of the 
energy dissipation field has such dependence: 

K ee (R) ~ e 2 [L/R)‘ 2< ’ 2 ~ C ’ 4 ‘(R/ 7y) 2 ( A_Ac ) (1.7) 

where A c = 2 — £ 2 . It has been argued |l4|] that if 
A < A c (a situation referred to as the “subcritical sce¬ 
nario”), then K41 is asymptotically exact for infinite Re. 
Then 2£ 2 = £4 = 4/3 and the outer scale L disappears 
from & In that case the exponent /i is identified with 
2(4/3 — A), and the renormalization length is the inner 
length rj. In fact, it will be shown here that the exponent 
A can be computed explicitly, and that it takes on ex¬ 
actly the value A = A c . Asa result of this the correlation 
K ee (R) can be shown to depend on R like 


where P is the transverse projection operator P = 
—A _2 V x V x . It was explained in Paper I that develop¬ 


ing a perturbative approach [|l2 Ul7||l8| . |l9|l for the correla¬ 


tion functions and response functions in terms of the Eu- 
lerian velocity u(r, f) results in a theory that is plagued 
with infra-red divergences. On the other hand one can 
transform to new variables, and after the transformation 
(which amounts to infinite partial resummations in the 
perturbation theory) one find a renormalized perturba¬ 
tion theory that is finite, without any divergences in any 
order of the expansion (cf. JnJ and I) . The new vari¬ 
ables are obtained from the Belinicher-L’vov transforma¬ 
tion Hell. 


r[r 0 \r,t\ = u[r + p(r„,t),t)] , 


( 2 . 2 ) 


where p(tq, t) is the Lagrangian trajectory of a fluid point 
started at point r = r 0 at time t = to 


K ee (R) =~ e\R/L) 


A4-2C2 


( 1 . 8 ) 


In other words, the critical situation A = A c results 
in the disappearance of the inner renormalization scale 
and the appearance of the outer renormalization scale in 
(1.8). In addition one notes that K ee (R) decays as a func¬ 
tion of R (i.e. the correlation is mixing) only if £4 < 2/2 
which implies deviations from K41. Thus we will argue 
in paper III that the critical scenario A = A c goes hand 
in hand with multiscaling if K ee (R) is mixing and then 
p, is identified with (4 — 2 / 2 . 

The quantity that displays the existence of the anoma¬ 
lous exponent A in the simplest possible way is the non¬ 
linear (4-point) Green’s function. Thus, after reviewing 
in Section 2 some essential results from previous work, we 
turn in Section 3 to the nonlinear Green’s function. We 
develop its diagrammatic representation and show how 
all the diagrams in its expansion can be resummed ex¬ 
actly into ladder diagrams. We discuss the properties of 
the ladder diagrams, and demonstrate in Section 4 that 
the dressing of the ladder diagrams does not change the 
fact that they resum to give power laws with anomalous 
exponents. In section 5 we firstly demonstrate that the 
resummed equation for the Nonlinear Green’s function 
indeed has a nonperturbative solution with an anoma¬ 
lous exponent A, and then derive a sum rule that allows 
us to conclude that A = A c = 2 — £ 2 - Last Section 6 
concludes this paper. 


2. SUMMARY OF PREVIOUS RESULTS 

The starting point of the analysis are the Navier-Stokes 
equations for the velocity field of an incompressible fluid 
with kinematic viscosity v which is forced by an external 
force f(r, t): 

u + P(u • V)u = Pf (2.1) 



P( r o,t)= u[r + p(r 0 ,r),r]dr . (2.3) 

Jo 

The natural variables for a divergence free theory are the 
velocity differences 

w(r 0 |r,f) = v[r 0 |r,t] — v[r 0 |r 0 ,t] . (2.4) 


Since the averages of quantities that depend on one time 
only can be computed at t = 0 , it is clear that the average 
moments of these variables are the structure functions of 
the Eulerian field: 


S n (\r- r 0 |) = (|w(r 0 |r,f)| n ) . (2.5) 

It was shown 0 that these variables satisfy the 
Navier Stokes equations, and that one can develop 
(cf. I) a perturbation theory of the diagrammatic 
type in which the natural quantities are the Green’s 
function G Q| g(ro|r, r', t, t') and the correlation function 
F a p(r 0 \r,r\t,t'): 


G a p( r 0 |r, r ' ,t,t’) 


F a p{ r 0 \r,r',t,t') 


<5(w a (r 0 |r, t)) 

SMr'ff’) h ^ Q 
{w a {r 0 \r,t)wp{r 0 \r',t’)) . 


( 2 . 6 ) 

(2.7) 


In stationary turbulence these quantities depend on t' — t 
only, and we can denote this time difference as t. The 
quantities satisfy the well known and exact Dyson and 
Wyld coupled equations. The Dyson equation reads 


— - i/AG Q/3 (r 0 |r,r',t) = G^rolr,r', 0 + )6(t) 

+ /dr 2 GO,(ro|r,r 2 ,0+) j dr 1 J dt± (2.8) 
x E5 7 (r 0 |r 2 , it, fi)G 7/ a(ro|ri,r', t — ti ), 


where G° /3 (r 0 |r,r', 0 + ) is the bare Green’s function de¬ 
termined by Eq.(3.20) in paper I. We will only need 
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the asymptot ic properties of this function, given in 
Eqs.( ^.17 2.1S ). The Wyld equation has the form 


F a p(r 0 \r, r',t) = j dvidr 2 J dt 1 dt 2 G ol y(r 0 \r,r 1 ,t 1 ) 
Dys (ri - r 2 ,t- t x +t 2 ) + <f> 7 5(r 0 |ri,r 2 ,t - G +t 2 ) 


x Gs(3{ro\r',r 2 ,t 2 ) . 


(2.9) 


In equation ( |2.8| ) the “mass operator” E is related to the 
“eddy viscosity” whereas in Eq. ( |2.9| ) the “mass operator” 
$ is the renormalized “nonlinear” noise which arises due 
to turbulent excitations. Both these quantities are de¬ 
pendent on the Green’s function and the correlator, and 
thus the equations are coupled. The diagrammatic nota¬ 
tions notations for the Green’s function, double velocity 
correlator and vertex are presented in Fig.p]. The dia¬ 
grammatic representation of E and $ is well kn own , and 
i s re produced in Fig.|j. In deriving equations (|2.8| ) and 
(2.9) one assumes that the driving force of the equation 


for w(ro|r, t) is a Gaussian random force whose covari¬ 
ance is the factor Dys which appears in (2.S). For future 
reference we need to remind the reader of the existence of 
a so called “principal cross section” in the diagrams for 
$. This is defined as a section of the diagram that divides 
left and right parts of the diagrams by cutting through 
wavy lines of double correlators only. Each diagram in 
the series for has a unique principal cross section, and 
see paper I for more details. 

The main result of I is a proof of the property of “lo¬ 
cality” in the Dyson and Wyld equations. This property 
means that given a value of |r — ro| in Eq. (|2.8| ), the im¬ 
portant contribution to the integral on the RHS comes 
from that region where (iq — r 0 | and |r 2 — r 0 | are of the 
order of |r — ro| . Moreover, the diagrams in the expan¬ 
sion for E, see Figj2[ are also local in the same sense: all 
the intermediate coordinates in all the integrals must be 
of the same order of magnitude, i.e. |r — ro |, in order to 
give an appreciable contribution. In other words, all the 
integrals converge both in the upper and the lower limits. 
The same is true for the Wyld equation, meaning that in 
the limit of large L and small 77 these length scales dis¬ 
appear from the theory, and there is no natural cutoff in 
the integrals in the perturbative theory. In this case one 
cannot form a dimensionless parameter like L/r or r/rj 
to carry dimensionless corrections to the K41 scaling ex¬ 
ponents. For 77 <C |r — r 0 1 <1 scale invariance prevails, 
and one finds precisely the K41 scaling exponents: 


G a/ 3 (Ar 0 |Ar, Ar',A z f) = A 52 G a/3 (r 0 |r, r', t), 
F a / 3 (Ar 0 |Ar, Ar',A 2 f) = A f2 F Q/3 (ro|r, r', t) . 


( 2 . 10 ) 


One can derive two scaling relations which hold order by 
order, i.e 


Note that nonperturbative effects may change these 
scaling relations. In this paper we will see that even if 
we use the K41 values of these scaling exponents, we dis¬ 
cover anomalous scaling in turbulence for higher order 
objects. We need to look for them using nonperturbative 
methods. The aim of this paper is to present a careful 
study of the nonlinear Green’s function in which anoma¬ 
lous exponents appear in a most transparent way. 

Some additional results that were proved in I and that 
we need for our analysis below have to do with the time 
integrals and time scales of the Green’s functions in the 
asymptotic regimes r>r' and r <C r'. We find that the 
time integral of the Green’s function satisfies 

dt G(0|r, r', t)) = r(r, r')G° (0|r, r', 0+) (2.12) 

where 

r(r, r') ~ r(r)(r , /r)“ for r r (2-13) 

with some scaling exponent a > 0 , and with 

r(r) ~ [e]" 1 / 3 r 2/3 . (2.14) 

In the standard phenomenology of turbulence the time 
scale r(r) is known as the characteristic turn-over time 
of an r-eddy. We showed that this time scale appears 
naturally in a scaling relation, i.e. 

r(r) = r/\JS 2 (r) . (2.15) 

In the opposite limit, i.e. for r <C r', we found that 

r(r, r') ~ r(r)(r'/7') /3 for r <C r' (2-16) 



where f3 < 0. The zero-t ime Green’s function 

G°(0|r, r', 0 + ) which appears in (2.12) and serves as an 


initial condition for the Dyson equation has the following 
asymptotic properties: 

G^(0|r,r',0 + )~(l/r ') 3 for r» r ', (2.17) 

Ga/?(0|r, r', 0 + ) ~ r-r'/r' 5 for r' > r . (2.18) 


Finally, we need to quote also the results for the asymp¬ 
totic properties of the correlator which were derived in I. 
For r « r' 


F af} { 0|r, r', 0) ~ e 2/3 (r 2 / 3 + r/r ,1/3 ) . (2.19) 

The correlator is symmetric in r and r' and the refor e in 
the opposite limit one just replaces r and r' in (2.19). 


3. THE NONLINEAR GREEN’S FUNCTION 


2 z + £2 — 2 , z -j- 2 C 2 — 2 . 


( 2 . 11 ) 


In this section we begin to expose a mechanism for 
anomalous scaling which operates in the context of 


The solution is z = ( 2 = 2/3. It was also shown that 
£2 = -3. 
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various many-point objects which appear in the the¬ 
ory. Such objects depend on two or more space coor¬ 
dinates. For pedagogical purposes it is convenient to 
discuss one of the simplest objects which display anoma¬ 
lous behavior, which is the nonlinear Green’s function 
(ro|a?i, x 2 , £ 3 , X 4 ) defined as, 


G% 0 ' lS (r o \x 1 ,X2,x 3 ,X4) 


£w a (r 0 |xi) fa> 7 (r 0 |x 2 ) \ 
5hp{r 0 \x 3 ) 6 hs(r 0 \x 4 ) / ’ 


(3.1) 


where for brevity we use the notation Xj = {r 
In a Gaussian theory (which our is not) this quantity 
would be the product of the linear Green’s functions, 
G a/ 3 (r 0 \xi,X 3 )G^ s (rQ\x 2 ,X 4 ). In a non-Gaussian theory 
it is natural to assume that this quantity is a homoge¬ 
neous function of its arguments when they are in the 
scaling regime. This means that 

G^'VolAri, A 2 G, Ar 2 , \ z t 2l Ar 3 , A 2 t 3 , Ar 4 , A 2 f 4 ) 

= \ u G% / 3 ' l 5 (r 0 \x 1 ,x 2 ,x 3 ,X 4 ) (3.2) 


From the Gaussian decomposition of this quantity we 
would guess that £4 = 2 £ 2 = — 6 . The proof of local¬ 
ity in I means that there is no perturbative mechanism 
to change this scaling index. On the other hand, this 
quantity, which is a function of four space-time coordi¬ 
nates Xi has scaling properties that are not exhausted 
by the overall scaling exponent £ 4 . We will show that 
when we consider its dependence on ratios of space-time 
coordinates in their asymptotic regimes we pick up a set 
of anomalous scaling exponents. Our first objective is to 
show that in the regime r± ~ r 2 <C r 3 ~ r 4 the diagram¬ 
matic expansion of this object produces logarithms like 
ln(r 3 /ri) to some power. Next we will prove that the 
sum of such logarithmically large contributions is given 
by (r 3 / 7 'i) A with some anomalous exponent A. To make 
the appearance of anomalous exponents evident we begin 
with the simplest object that resums to logarithms, i.e. 
the series of “ladder diagrams”. 


A. The ladder diagrams 

The diagram matic representation of the nonlinear 
Green’s function (0 is obtained as follows. In the spirit 
of the Wyld expansion 0 © one can begin with the di¬ 
agrams for w which are shown in Fig.|3|. After differen¬ 
tiating with respect to a force we get the diagrams for 
the unaveraged response 8 w a (rQ\x\)/5hp(r^\x 3 ) shown 
in Fig§ We recall (see I) that each of these diagrams 
has a principal path of Green’s functions connecting an 
entry denoted by a wavy line to an exit denoted by a 
straight line. At this point we can take any combination 
of two diagrams, one for 5w a {ro\xi)/5hp(ro\x 3 ) and one 
for 6 w a (ro\x 2 )/Shp(ro\x 4 ,), glue them together according 
to the Gaussian rules and then perform the Dyson-Wyld 


line resummation. Clearly every resulting diagram has 
two principal paths, one going from xi to x 3 and the 
other from x 2 to X4, see FigJ|a. Every diagram begins 
with two wavy lines and ends with two straight lines. 
One infinite sum of contributions results from having all 
the averaging done in each tree separately. Such a sum 
results in the first diagram on the RHS of Fig. 4a, which 
is precisely the product of two dressed Green’s functions. 
Next we have infinite number of diagrams in which these 
two principal paths are connected via fragments that re¬ 
sum to a dressed correlator. This is the sum that gives 
rise to the second diagram on the RHS of Fig. 4a. Fol¬ 
lowing are infinite sums that result in connection via two, 
three etc. dressed correlators. The sum of terms appear¬ 
ing in Fig||a is referred to as the sum of “simple ladder 
diagrams”. We note that this sum can be represented ex¬ 
actly in graphical notation as shown in Fig.|5p. We used 
here explicitly the simple topology of the simple ladder 
diagrams. 

In the simpler case of the passive scalar that was pre¬ 
sented in [Q the simple ladder diagram tell the whole 
story. This is not the case here. There exist an infinite 
class of diagrams that decorate the simple ladder dia¬ 
grams as shown in Fig.|(i|. The result of the summation of 
all the simple ladder diagrams in Fig.|| which is denoted 
as a thin circle is added to the diagrams that contain 
decorations of the vertices, see Fig H These decorations 
are identified as contributions that can be resummed to 
dress the vertices. In these diagrams there are two types 
of dressed vertices that we denote as A and B respec¬ 
tively. The decorated vertex A has one straight tail and 
two wiggly tails like the bare vertex, and therefore A can 
be considered as the dressing of the bare vertex. The ver¬ 
tex B has two straight tails and one wiggly, and it does 
not have a bare counterpart. 

In the general theory one finds an infinite series of di¬ 
agrams that dress the vertices A and B, and this series 
is shown in Fig.|?]. It is therefore clear that the sum of 
all the diagrams appearing in Fig.|5| and in Fig.^j can be 
represented as diagrams that have the topology of the 
simple ladder diagrams but with dressed vertices. These 
diagrams are shown in FigJ|. We refer to the diagrams 
shown in panel a of Fig|| as the sum of “dressed simple 
ladder diagrams”, and denote the result of this summa¬ 
tion with a bold circle. We note that one cannot have 
two vertices of type B in one rung of the ladder. Such 
a rung will call for a connection via a long straight line, 
an object that does not exist in this theory. Panel b of 
Fig.|| represents the exact resummation of all the dia¬ 
grams having the topology of the dressed simple ladder. 
Note that the resummed series which is represented as a 
heavy circle appears on both the LHS and RHS of this 
equation. This is not yet the fully resummed nonlinear 
Green’s function since there are additional diagrams with 
different topology. 

Next we need to discuss the diagrams whose topology 
differs from that of the dressed simple ladder diagrams. 
Such diagrams are obtained when the connections be- 
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tween the principal paths intersect, or when there appear 
decorations that cannot be absorbed into the dressed ver¬ 
tices. The simplest examples are shown in Fig. |[i, in 
which the vertices have been already dressed by summing 
the corresponding decorations of the diagrams shown. 

The resummation of all the diagrams for the Nonlin¬ 
ear Green’s function is shown in Fig.[|b, where the series 
for the renormalized rung is given in Fig. 9a. The di¬ 
agrams in the series can be classified into two classes, 
which we refer to as “two-eddy reducible” and “two eddy 
irreducible” diagrams. To define these we need first to 
define “interior” Green’s functions as Green’s functions 
belonging to one of the two principal paths but not to an 
entry or an exit. “Two-eddy reducible” diagrams are dia¬ 
grams that can be cut into two pieces by cutting through 
two interior Green’s functions belonging to two principal 
paths without cutting any other line. “Two-eddy irre¬ 
ducible” diagrams are diagrams that are not two-eddy 
reducible. As examples consider the diagrams in Figj5|. 
In 4a the two first diagrams are two-eddy irreducible, and 
the next two are two-eddy reducible. In Fig.^j diagram 
1 is two-eddy reducible whereas diagram 2 is two-eddy 
irreducible. 

The sum of all the two-eddy irreducible dia¬ 
grams gives rise to the “two eddy mass operator” 
S 2 (ro|aq, x 2 , X 3 , X 4 ). This object is given in terms of the 
infinite expansion shown in FigJlC|a. Using this object 
we can sum all the diagrams belonging to the nonlinear 
Green’s function, as shown in Figjujb. This is an ex¬ 
act equation for the nonlinear Green’s function, which 
is the direct analog of the Dyson equation for the usual 
Green’s function. Note that the fully resummed non¬ 
linear Green’s function is denoted by a gray circle with 
bold circumference. Thus in Fig. 8a the bold empty 
circle which is the sum of all the dressed simple ladder 
diagrams is added to all the diagrams that have more 
complex topology, to yield the gray circle object. Note 
that the equation in Fig.|Ic|) reduces to the equation in 
Fig.^b if we substitute the first three terms of Fig -0a 
instead of the two eddy mass operator in Fig.|I(jjb. 


B. Logarithmic contributions in simple ladder 
diagrams 


In this sub-section we demonstrate that the leading 
contribution to the ladder diagrams involves logarithmi¬ 
cally large terms 12 TTi||n| . [T^ . [i7i| . po| . To this aim we in¬ 
troduce the time integral of the nonlinear Green’s func¬ 
tion defined by 


linear Gre en’s functions and is denoted as diagram 1. 
Using Eq.( [2.12 ) we compute 


G 2*,o( r oki,a ; 2,r 3 ,r4) (3.4) 

= r(r 1 , r 3 )r(r 2 , r 4 )G 0 (r 0 |r 1 , r 3 ,0 + )G o (r o |r 2 , r 4 ,0 + ) . 


Con sider now the regime rq 
Eq.([2.13) to estimate 


r" 2 <C rq ~ r 4 . We can use 


G 2*o( r o|zi 5 Z2,r 3 ,r 4 ) 


T(ri)r(r 2 ) 

(r 3 r 4 ) 3 



Consider now the next term which is the shortest lad¬ 
der in Fig. |B| (i.e the ladder with one rung and two ver¬ 
tices which is labeled as diagram 2). It is denoted as 
G^ 2 (r 0 |oq,a; 2 ,r 3 ,r 4 ) and it is written as 


G 2 * 2 ( r oki, ® 2 , r 3 , r 4 )~ J dXidxjG(r 0 \xi,Xj)G(r 0 \x2,Xi) 

-F(r 0 \xi,Xj) J dt 3 dt4G(r 0 \xj,x 3 )G(r 0 \xi,X4) . 

(3.6) 


d__d_ 

dri dr j 


In writing Eq.(3.6) we dropped the tensor indices from 
G and F. The range of integration that is of interest for 
us is the range rq ~ r 2 -C ry, ry <C r 3 ~ r 4 . In t his r ange 
we can use Eq.(2.13) and evaluate the integral (3.6) as 


G 2* 2 ( r o|zi,a:2,r 3 


^ r(n)r(r 2 ) / rqr 2 \ P + 1 
1-4 (rqr 4 ) 3 Vr 3 r 4 / 


dr 


dr 


— / ^F(r 0 |ry,ry)T(ry)T(ry) 


r 


/ r 2 

7 V2 ' J 


(3.7) 


(provided rq ~ r 2 <C ry, rj «r 3 ~ r 4 ) . 


To proceed we will use the K41 scaling exponents to 
evaluate the correlators. In other words we evaluate 
F(r 0 |ry,ry) as min{ry 2 / 3 , ry 2 / 3 } according to Eq.( |2.19 ). 
This is not an exact step, and is used to demonstrate in 
the swiftest way the existence of logarithmic divergences 
which lead to anomalous scaling. In Section 6 we will 
derive an exact sum rule that will establish the existence 
of the anomalous exponent and will determine its value. 

Consider now separately the two possibilities ry < rj 
and ry > rj. In the former case we find that the ry 
integral contributes mostly in the upper limit and the rj 
integral in the lower limit. In the latter case the situation 
is inverted. In both cases the integrals have their largest 
contribution in the regime ry ~ ry. We thus can evaluate 
© as 


G 2 *(ro| 2 q,£ 2 ,r 3 ,r 4 ) = / dt 3 dt 4 G 2 (r 0 |aq,z 2 ,x 3 ,z 4 ) . 

(3.3) 


G 2 t , 2 ( r o|2 ; i,Z2,r3,r 4 ) 


G “o( r o|zi,Z2,r 3 , 



dry 

ry 


= A 0 G 2 ‘ 0 (ro|aq, a; 2 , r 3 , r 4 ) ln(r 3 /rq), (3.8) 


The zero order contribution to this quantity is obtained where Ao is a dimensionless constant of 0(1). The next 
from the first term in Fig||i which is the product of two order contribution is 
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G^roki, x 2 ,r 3 , r ±)~JdxidxjG(r 0 \xi,Xj)G(r 0 \x 2 ,Xi) 
d d f 

X ~Q^r Qj. n**\*U ^0 ) J dXndXmG (y q \ Xj , X n ^G(Y 0 \ X{ , X m ) 

d d f 

x——-— -F(ro|a; n , x m ) dt 3 dt4G(jo\x n , x 3 )G(ro\x m , X4). 


dr n dr m 


(3.9) 

It can be seen in analogy to the sit uation in f3) that 
the intermediate integrations in (3.E) peak in the regime 
Ti ~ Tj and r n ~ r m . The largest contribution to the 
integral then comes from the regime rq ~ r 2 -C rq ~ 
rq <C r n ~ r m <C r 3 ~ r 4 - In this regime the integral is 
evaluated as 


G “4( r oki.a2,r3,r 4 ) 

r 

G 2 *o( r ol a; i>* 2 )r 3 ,r 4 ) / 

J n 

G 2 t .o( r o|a:i ! *2,r3,r4)- A 0 In 


dri 


Ti 


dr n 
r n 

1 2 


(3.10) 


Note that we have asserted here without proof that the 
coefficient in fron t of the logarithm is the square of the 
coefficient in (3.8). This is intuitively acceptable because 
of the repetitive nature of the ladder structure. It is 
however an important point and therefore we will prove 
it in section 5 by analyzing the resummed equation for 
the nonlinear Green’s function. 

It becomes believable now that the ladder with n rungs 
contains a contribution of order [A 0 ln(r 3 /ri)] /n\. We 
reiterate that there are additional contributions from 
other regimes of the intermediate integrations and they 
will contribute lower powers of logarithms. One can take 
them into account, but this only serves to renormalize 
the value of Ao, as is proven below. At this level of 
demonstration it is enough to take the leading order con¬ 
tribution from each n-rung ladder and notice that the 
summation of all these contributions will give a term pro¬ 
portional to (7'3 /Vi) A °. 

Next we need to prove that the decorations of the sim¬ 
ple ladder diagrams do not change the qualitative picture 
of anomalous scaling discussed above. 


4. PROOF OF RIGIDITY 

A useful concept in the demonstration of anomalous 
scaling is the concept of rigidity, which is an order by 
order property of diagrams that we are going to use re¬ 
peatedly. We explain the concept with the help of Fig.pd|. 
Consider the diagram in Figjll|a. The fragment within 
the dashed circle is a part of the mass operator E. Sup¬ 
pose that x a is smaller than Xb■ the property of locality 
which was proven in I means that all the significant con¬ 
tributions to the diagrams come from the range of inte¬ 
gration in which x 3 ,x 2 ,x 3 and X 4 are all lying between 


x a and Xb- The property of rigidity is stronger, and it re¬ 
lates to fragments of diagrams whose exterior points are 
vertices. For the diagram in Fig.|ir|a rigidity says that if 
we fix the coordinate xi then for x a -C X\ and Xb x\ 
the largest contribution to the diagram comes from the 
regime of integration x 2 ~ x 3 ~ X 4 and all these space- 
time coordinates are of the order of X \. We call this prop¬ 
erty “rigidity” to give the intuitive feeling that one can 
stretch at will the diagram such that \x a —X\ | and \xb~x\ \ 
become very large, but because x\ was determined, the 
positions x 2 ,x 3 , X 4 are rigidly fixed to the vicinity of X\. 
We call the diagram rigid if it has this property with re¬ 
gards to fixing either x\ or x 2 . In diagrams that have 
three external legs to them, like the one in Fig. 0b which 
belongs to the series of 3-point correlator, rigidity of the 
diagram implies that the inner fragment is rigid with re¬ 
spect to fixing either X\ or x 2 , or x 3 . Thus if we fix x 3 
than x a ,Xb and x c can be brought arbitrarily far away 
from x\, and yet the main contribution to the integral 
comes from the regime x 2 ,... ,x-j are all of the order of 
xi, etc. This definition of rigidity extends to diagrams 
with an arbitrary number of external coordinates, like in 
Fig. 10c. Such a diagram is called rigid if fixing any 
of the coordinates x 3 ,...,xe results in the main contri¬ 
bution to the diagram coming from the regime in which 
the rest of these coordinates are of the same order as the 
fixed one, independently of the positions of the external 
coordinates x a , ■ ■ ■, Xf. 

Lastly we want to clarify the concept of a “stretched” 
diagram. A stretched diagram is a a diagram in which 
there is a definite ordering in the positions of the vertices 
that are being integrated upon. Consider for example the 
diagram in Figjlljc. Suppose that the positions r 2 , r 3 , r 4 
of the group of vertices x 2 , x 3 ,X 4 are larger than the 
positions ri,re,r 3 of the group xi,xe,x 3 . We call this 
diagram stretched if all the positions r Q ,r/ and r e are 
smaller than any of ri,r 6 ,r 3 , and the positions r&, r c , rd 
are all larger than r 2 ,r c ,rd- Rigidity of a fragment will be 
intuitively understood as the resistance to stretching. In 
other words, one can think of the propagators as springs 
that are being pulled to stretch the diagram. If it is suf¬ 
ficient to fix one coordinate of the fragment, in this case 
any of the coordinates X\ — xq such that the main con¬ 
tribution to the diagram comes from the regime that all 
the coordinates in this group are of the same order, we 
call the diagram rigid. 

We prove now that the ladder diagram are rigid. 


A. Rigidity of the dressed vertices 

The aim of this sub-section is to answer the following 
question: how should we evaluate a dressed vertex in a 
stretched diagram. Every dressed vertex has three coor¬ 
dinate designations, see for example Fig.^j, in which these 
coordinates are x a ,Xb and x c respectively. For fixed co¬ 
ordinates Xi to X 4 and Xb,Xd we need to integrate over 
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a : a , and x c . We will show here that the main contribution 
to the integral comes from the region x a ~ Xb ~ x c . As 
a result we will be able to estimate the dressed vertex 
as 1/rt,. Finally this estimate will allow us to repeat the 
argument that rf, contributes mostly in the vicinity of rd- 
Thus the estimate leading to the logarithmic divergence 
will remain unchanged. 

We begin by considering the various possible integrals 
that depend on the local geometry around a given vertex 
in a stretched diagram. These are the local geometries 
shown in Figs.|l2] and [l^. In these figures the dashed 
line represents either a straight or a wiggly half line. We 
think of these fragments as parts of a bigger diagram in 
a “stretched” configuration, meaning that in the larger 
diagram there are no turn backs in the positions of other 
vertices that are not shown in Figs|l2] and |u| In other 
words, if we consider Figjijja, the positions of all the ver¬ 
tices attached to x\ are smaller than n, and the vertices 
attached to x 2 and 23 are all larger than r 2 and 23 . For 
lib it means that all the vertices attached to x\ and 23 
have positions smaller than j*i and 23 , and all the vertices 
attached to 2 2 have positions larger than r 2 . 

The question that we want to answer is the following: 
when we integrate over f, f 2 and t^ where is the dominant 
contribution to the integral over r. In order to answer 
this question efficiently we decompose the various inte¬ 
grals to their elements, which are shown in Fig.[n|. Every 
vertex is a junction of two wavy lines and one straight 
line. The Green’s function belonging to the latter is con¬ 
sidered together with the vertex as one element, shown 
in panels a and b. The wavy line can belong to either a 
correlator considered in panels c and d, or to a Green’s 
function considered in panels e and f. In panels c-f the 
vertex is excluded from considerations, as is indicated 
by the dashe d lin e. We eval uate the time integrals on 
the basis of ( [2.13 ) and (2.16). For brevity we will take 
a = (3 = 0. This means that our results shown in Fig.13 
are bounds, and the actual situation is always better in 
terms of the conclusion where the integral over r con¬ 
tributes most. 

The estimate of Fig.^dja is r~ 5 which is obtained from 
the estimate of the vertex like 1/r and of the Green’s 
function like ri/r 4 . The integral over t is restricted by 
r(ri) and does not contribute. Fig. & is estimated as 
r~ 4 which stems from the estimate of the Green’s func¬ 
tion as 1/r 3 . Fig. Ef is r° because the correlator is 2 j 2 / 3 
and the time restriction is r(rj ). On the contrary Fig.[14|d 
is r 4 / 3 because both the correlator and the time restric¬ 
tion contribute r 2//3 each. Figjl4[; is again r° because the 
Green’s function and the time restriction are determined 
by rj . On the contrary 13f gives r 5 / 3 due to r/r ^ 4 in the 
Green’s functions and r 2 / 3 from the time restriction. 

To evaluate the integral over 2 in Fig JJa we need to 
combine Figjl^a with Fig.[i~i)d and Figjl4f in all possible 
combinations, upon performing the integral over r we 
see that we always have the contribution coming form 
the upper limit,i.e. when r is of the order of r 2 or 23 . 


or 


To evaluate Fig.|l2|b we combine Figjl4|a with FTgJTdj: 
Kig.[T~i|e and with Fig.[TT|d or Fig.[T~i[f. Here the resulting r 
integral contributes mostly in the lower limit, when r is 
or the order of rq or 23 . Similarly we find that in 11c the 
major contribution is in the lower limit, as is the case for 
Fig.12c. In Fig.|l3|a and Fig.|l3|b the major contribution 
is in the upper limit. This information is summarized in 
Fig.[l5| that shows where the integral over 2 contributes 
for any possible arrangement of propagators. The arrow 
shows whether the major contribution is coming from 
the upper or the lower limit. We again want to invoke 
the mechanical analogy of springs, and point out to the 
reader that one can guess where the major contribution 
comes from simply by counting how many springs pull 
in each direction. We note that the “spring strength” is 
not the same for a correlator or a Green’s function. But 
in any combination two springs are always stronger than 
any single spring. 

These results can be used now to prove immediately 
that the simply decorated vertices in Fig|(] have their 
main contribution (in the stretched situation) in the re¬ 
gion of integration x a ~ Xb ~ x c . Consider for example 
the decorated vertex A in the diagram numbered as (1). 
Let us stretch the diagram such that 21 and 2 2 are the 
smallest positions and 23 and 24 the largest. Fix 2 ;, and 
realize that the results summarized in Fig.^J imply that 
the main contribution is when r a is of the order of 2 ;,. 
Similarly, the integral over r c will contribute in its lower 
limit, again for r c or the order of 2 (,. Using no w the 
proven locality x a ~ Xb ~ x c , the scaling relations ( 2 . 11 ) 
immediately imply that the evaluation of the decorated 
vertex cannot be changed from the evaluation of the bare 
one, and it is therefore 1 / 2 ;,. 

We could go on to any chosen decoration of the ver¬ 
tex and reach a similar conclusion. Consider for example 
the decoration of the same vertex A shown in Fig.|ljJ) . 
Fix the position 23 for example. Stretch the springs con¬ 
nected to 2 q and 2 b such that r a < r\ < 23 < 2 2 < 2b. 
The spring analogy which says that two springs are 
stronger than a single one implies that the dominant con¬ 
tribution comes from 21 of the order of the smaller of 24 
and 26 - Similarly, 2 2 will be of the order of the larger of 
24 and 25. Next we repeat the argument at 24 or 2 6 and 
24 or 25 and conclude that they contribute in the vicin¬ 
ity of 23 . When all the coordinates are of the same order 
the scaling relations protect the evaluation of the dressed 
vertex and keeps it as I/23. The spring analogy helps to 
understand that when we decorate the vertex further the 
rigidity is even more pronounced since there are more 
springs to keep it in place. 


B. Rigidity of the Rungs 

Having demonstrated the rigidity of the 3-point ver¬ 
tices, we can proceed now to examine the rigidity of the 
4-point rungs. For the case of the dressed simple lad- 
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ders, the argument is immediate. Consider the diagram 
in FigJ^a, and focus for example on the diagram desig¬ 
nated as 1. Suppose that we stretched n and r 2 to be 
much smaller than r a and and r% and to be much 
larger than r a and tv Can we have a significant contri¬ 
bution from the region in which r a becomes much smaller 
than rb or vice versa? The analysis summarized in Fig.|l5| 
or the spring analogy imply that the answer in negative. 
Both vertices will contribute in the same neighborhood 
r a ~ rb- This is the property of rigidity of the dressed 
simple rung. 

Next we are going to consider the two-eddy mass oper¬ 
ator £ 2 ( 1*0 hi, X 2 , X 3 , £ 4 ) which is positioned in a ladder 
like in Fig.pp (diagram 3) with another such mass opera¬ 
tor above and below. Here “above” means that the low¬ 
est two coordinates of the next mass operator are much 
larger than x c , Xd, and “below” means that the upper 
two coordinates of the previous mass operator are much 
smaller than x a , Xb■ Consider this object with one fixed 
coordinate, say r a , and integrate over the other three co¬ 
ordinates. We will prove now that the main contribution 
to this integral comes from the regime in which all the 
three coordinates are of the order of the fixed one. 

To this aim consider Fig.p^a. We display here a block 
representing a dressed rung of the ladder all of whose 
coordinates are of the same order of magnitude, except 
for the two blobs whose coordinates are much larger or 
much smaller than all the other coordinates. These two 
blobs must be connected to the rung with at least two 
legs, else the diagram would be one-eddy reducible. 

Consider the case that the blob includes just one ver¬ 
tex, say £ 3 , and it is indeed connected via exactly two 
legs, say to the vertices X 5 and Xq. In this case we are in 
the situation of FigfTsfb, and the vertex will contribute 
mostly when X 3 is of the order max(a :5 , Xq). If the vertex 
gets decorated, the proof of rigidity of the 3-point objects 
guarantees that the dominant contribution to the integral 
comes from the regime in which the positions of all the 
vertices are of the same order, say 7-3 ~ rg ~ rio- The 
scaling relations imply that this vertex can be evaluated 
as, say, l/ra, exactly like the bare one. Consequently the 
blob contributes mostly in the vicinity of max(a ;5 , x§) 
again. 

Next we discuss the effect of adding additional legs to 
the connection of the blob to the main block of the rung. 
These legs can be either correlators F(r a \x u ,Xd) or one 
of the Green’s functions G(r 0 \x u ,Xd) and G(r 0 \xd, x u ). 
Here x u belongs to the upper block x u ~ £ 3 , and Xd 
belongs to the lower block Xd ~ X5. We will discuss the 
effect of the addition of one such leg, and show that it 
makes the integration over X3 peak more strongly at its 
lower boundary. The addition of more legs just enhances 
this tendency. 

Adding a correlator F(r n \x u ,Xd) requires the addition 
of two vertices at x u and Xd, two Green’s functions, and 
two space time integrations over x u and Xd- This leads to 
an additional dimensionless factor (rd/r„) 7 . To evaluate 
7 note that upstairs we have r u 3 from the space inte¬ 


gral, 1 / r u 3 due to the Green’s function, and l/r u due to 
the vertex. The time integral upstairs contributes actu¬ 
ally r(rd) due to the time restriction on the correlator. 
Downstairs we have rd 3+2 ^ 3 from the space time integral, 
1 /rd 4 from Green’s function and vertex, and r ^ 2 / 3 as an 
evaluation of the added correlator. Consequently we find 
in this case 7 = 1 . 

If the added leg is G(r n \x u ,Xd) the situation up¬ 
stairs changes since the time integral is now bounded 
by r(r u ,rd ) < r(r u ). The situation downstairs also 
changes. We need to insert an additional correlator (in¬ 
stead of Green’s function) which contributes Td 2 ! 3 . The 
G{ro\x u ,Xd) contributes l/r^ 3 , and the space time inte¬ 
gral rd 3+2 ^ 3 . By power counting we conclude that in this 
case 7 > 1/3. 

Finally the added leg can be G(r 0 \xd,x u ), which is 
evaluated as rd/r u A . In this case we add upstairs a cor¬ 
relator which contributes r u 2 ^ 3 , and the space integral 
over the vertex gives r u 3 . The time integral upstairs in 
bounded by -r(rd). Downstairs we add a Green’s function 
evaluated as 1 /rd 3 , and the space time integration over 
the vertex contributes r^ +2 ^ 3 . In summary we compute 
7 > 4/3. 

We see that in all cases 7 is positive, leading to a 
stronger preference of the integral for peaking downstairs, 
as claimed above. This result is again in conformity with 
the intuitive picture that considers the propagators as 
springs that pull fragments of a diagram to each other. 
Adding more propagators to the detached blob makes the 
situation even more rigid. 

The analysis of the blob with smaller coordinates par¬ 
allels the above analysis exactly. Again the blob must be 
connected to the main rung with at least two propaga¬ 
tors, and the spring analogy works perfectly. 

Finally consider the situation in Fig.fl 6 |b, in which two 
blocks involving arbitrary resummations of two-eddy ir¬ 
reducible diagrams which are separated by at least three 
legs as shown. The two outer Green’s function are 
parts of the principal paths, and the dashed line repre¬ 
sents at least one additional connection, which is either 
F(jo\x u , Xd) or one of the Green’s functions G(ro|a; u , Xd) 
and G(ro|.Td, x u ). If all the coordinates were of the same 
order, this diagram would represent one rung of the lad¬ 
der, and its evaluation would be, according to the scaling 
relations, just like the evaluation of a bare rung. The aim 
of the present discussion is to show that indeed in the 
stretched diagram the largest contribution comes from 
the regime in which all the coordinates are of the same 
order. To show this we will assume that the two blocks 
have widely separated coordinates. All the upper coordi¬ 
nates are of the order of x u . and all the lower coordinates 
are of the order of Xd ■ We can fix one coordinate at will, 
say X\. The claim is that the major contribution to the 
diagram in Fig.|l 6 |b comes from the regime in which all 
the coordinates are of the order of x\. 

For pedagogical purposes it is useful to discuss first the 
two-eddy reducible topology in which there is no connec- 


tion (in addition to the two Green’s functions belonging 
to the two principal paths) between the two sub-blocks. 
If we assert that the coordinates in each block separately 
are of the same order, x u and Xd respectively, the integra¬ 
tion over all coordinates will result in a term proportional 
to j ri dr u /r u . Analytically this estimate is understood 
from the discussion of the bare ladder in Sect.^lr and from 
the rigidity of the dressed rung. In the spring analogy we 
have two springs below and two springs above, and this 
results in a neutral situation which is a “mechanical bal¬ 
ance”. Analytically this means that the integral is log¬ 
arithmic, having no preference to either upper or lower 
bound. Adding any one leg that turns these two sub¬ 
blocks into a two-eddy irreducible topology immediately 
introduces into this integral a factor ( rd/r u )' y a with 7 
taking one of the positive values found above. This factor 
forces the integral to peak in the regime r u ~ rd- Any 
additional leg connecting the two sub block only serves 
to enhance this tendency. The additional spring turns 
the whole rung rigid. 

One can come up with arbitrarily involved topologies of 
sub-blocks, say by introducing additional blobs at various 
coordinates between the two blobs shown in Fig fl6|b. The 
analysis presented here suffices to show that all these 
cases lead to the conclusion that the major contribution 
comes from the shell in space in which all the coordinates 
are of the same order. 


C. Dangerous contributions: ladders in ladders. 

The procedure outlined above consists of a proof of 
rigidity in every order of the diagrammatic representa¬ 
tion for the two-eddy irreducible mass operator which 
appears as a rung in the ladder. However, the anoma¬ 
lous exponents discussed above appears in a power law 
that results from the resummation of ladder diagrams to 
all orders. Therefore we need to be doubly careful about 
the appearance of resummations to all order within the 
blocks that represent the rungs in our ladder. In other 
words, we need to examine the possible appearance of 
ladders within the rungs of the ladders. 

In FigjT^ we present examples of two dangerous ladders 
in the rung of the ladder. In panel (a) the inner ladder has 
Green’s functions oriented in the same sense as the main 
ladder. In panel (b) the Green’s function are oriented in 
the opposite sense. In Figs.|l8|l7a and p~ 8 | b we present the 
inner ladders in r, t coordinates, to stress that the time 
restrictions on the Green’s functions force the two inner 
ladders to be oriented in the opposite direction in time. 

The understanding of the situation of Figjl7]a is not 
difficult. We have learned before that the insertion of a 
leg of G(r 0 \xd, x u ) results in a reduction factor of the or¬ 
der of (rd/r u ) 4 / 3 . The inner ladder in Fig.16a starts with 
two such legs, and this results in a total reduction of the 
order of ( Xd/x u ) 8//3 . On the other hand, the ladder itself 
will contribute an additional factor of ( r u /rd) A ■ We thus 


expect that in toto we will have a factor of (xd/x u ) 8 ^ 3 ~ A . 

This calculation show that if A exceeds 8/3 we will 
lose the property of rigidity and our theory will be in 
real danger. In fact we will show momentarily that the 
borderline of applicability of the theory is actually A = 
4/3. As long as we assume that the numerical value of 
A is smaller than 4/3 we can conclude that the inner 
ladder still acts as a spring which is at least as strong as 
an additional G(ro|x<i, x u ) leg. It is in no way dangerous. 

The understanding of the role of the second type of lad¬ 
der, Fig.|l7|b and |l 8 |b, calls for a different type of consid¬ 
eration. Suppose that we want to insert the ladder of 17b 
between two rungs (say, between the lowest and the next 
lowest rungs) of an outer ladder of the type 17a. How¬ 
ever 17b is oriented (by causality) upward to the right. 
Since the ladder in the 17a is oriented (again by causal¬ 
ity) upward to the left, inserting 17b between two rungs 
severely constrains the time range of the inserted ladder. 
This time constraint prevents a logarithmic divergence. 
In other words we cannot have a large factor of ( r u /rd) to 
some exponent in the inserted ladder. The lesson is that 
consequences of causality on the time restriction of the 
Green’s functions do not allow an insertion of one type 
of ladder in the other. 

However these are not the only dangerous situations 
that need further careful discussion. For example in 
Fig.^ we display yet another configuration of contribut¬ 
ing diagrams in which the bare Green’s functions that 
define the ladder are replaced themselves by ladders. In 
this case the estimate of the case shown in Fig. 16a is 
corrected by another anomalous factor of ( x u /xd) A ■ One 
can still state that for A < 4/3 the power —8/3 keeps 
the rung rigid, but we already come closer to a possible 
breakdown of the property of rigidity which would occur 
if A = 4/3. Indeed, we are going to show below that A 
attains precisely its critical value 2-/2 which for K41 
scaling is exactly 4/3. This fact means that nonpertur- 
bative effects which appear in infinite resummations are 
very important, and they may indeed lead to nontrivial 
renormalizations of the scaling exponents. 


D. Radius of the ball of locality 


Before continuing our study of anomalous scaling we 
pause to make use of the concept of rigidity to improve 
our estimate of the radius of the ball locality. The 
question arises in the context of the Green’s function 
G(0|r, r', t) and the correlation function F(0|r, r', t) when 
r and r' are of different order of magnitude. The result 
of this consideration will be that the radius of the ball of 
locality is always determined by the smaller of r and r'. 
The reader who is not interested in this issue is invited 
to go directly to Sect.|j. 

In the context of the Green’s function the question 
is what is the spatial domain which contributes mostly 
to the integration over ri and r 2 in Eq.(2.8). Due to 
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rigidity, in the stretched situation when r <C r' or vice 
versa the coordinates r\ and r 2 must be of the same or¬ 
der to contribute significantly. In this regime the mass 
operator E can be evaluated as (cf. paper I) S 2 {ri)/r\ 
~ e 2 / 3 r)~ 13 ^ 3 . The integrals over ri, r 2 and t\ are 
bounded f rfr(n)dri. Next the product of the two re¬ 
maining Green’s function can be always estimated as 
G(0|r, r', 0 + )/rf. The total T\ dependence of the inte¬ 
grand is r ^ 4 / 3 and therefore the integral contributes at 
its lower limit which is the smaller of r and r'. 

In the context of the correlator we need to analyze 
the integrals over ri, r 2 , t\ and t 2 in Eq. (2.9). In the 
stretched situation rigidity allows us to restrict our con¬ 
siderations to the regime n ~ r 2 . In this regime the 
mass operator $ can be evaluated (cf. I) as <S , |(n)/r 2 . 
The space-time integrals contribute to the ri depen¬ 
dence as before, because one of the time restrictions is 
r(min{r, r / }). Now the two Green’s functions are differ¬ 
ently oriented and can be evaluated as min{r, r'}/r j. I n 
toto we find the same evaluation of the integrand in ( p.9[ ), 
namely r 1 4//3 . Again the integral contributes at its lower 
limit which is the smaller of r and r'. 

The c onclusion is that due to rigidity, in both Eqs. 
( | 2 . 8 |j 2 .g| ) the relevant domain of integrations over ri and 
r 2 is n ~ r 2 ~ min{r, r '}. Now we can apply the prop¬ 
erty of locality proven in I to see that upon expanding the 
mass operators in their infinite series, all the intermediate 
coordinates have to belong to the same domain. Conse¬ 
quently, the global ball of locality has a radius min{r, r'}. 


5. ANOMALOUS SCALING 

In this section we discuss the resummed equation for 
the nonlinear Green’s function, and solve for the anoma¬ 
lous exponent that is associated with this function. The 
equation is shown graphically in Fig. [To)), and in analytic 
form it reads 


much larger than the inhomogeneous solution. In subsec¬ 
tion B we will evaluate the anomalous exponent exactly. 

In order to define the anomalous exponent A precisely 
we introduce the following function: 

T af 3 (ri,r 2 ,i?) = /^ 4 ^^r 0 |*i,* 2 , ! r 8 ,* 4 )< ) (* 3 , aS4 ) > 

(5.2) 

where (x 3 , X 4 ) is some function with a characteristic 
length scale R and characteristic time scale t(R). The 
anomalous exponent is defined via the limit 

lirn Vi • V 2 T’ Q , a (ri,r 2 , A) oc -i- . (5.3) 

V<ri,r 2 <CR Ti2 

The reasons for this somewhat cumbersome definition of 
A will become clearer in paper III of this series. 


A. Homogeneous Solutions: Qualitative Analysis 


The equation 0 ) is not an easy equation to solve. 
The nonlinear Green’s function is a function of four 
space-time variables, and it is a fourth-rank tensor. In 
order to gain insight on G 2 we will reduce it to a function 
of a smaller number of variables. Also, we are interested 
here in its scaling properties only, and we can simplify the 
discussion by dropping the tensor indices. This will allow 
us to develop a qualitative analysis that will dem onstrate 
the power law behaviour that is i mpl ied in (5.3). The ac¬ 
tual computation of the limit in ( p.3| ) will be done in the 
next subsection. 

The reduction in number of variables is done as follows: 
we integrate G 2 over the last two time variables t 3 and 
f 4 , and consider its value at t\ = t 2 = 0 and iq = r 2 = r, 
and r 3 = r 4 = R. We choose R 8 S> r. Define a new 
quantity < 72(15 R) via the equation 


^ 7<5 (ro I £ 1 , £ 2 , 2 ^ 3 , 24 ) = G“ 7 (ro|aq, a^G^rola^, x 4 ) 

+ f dx a dx b dx c dx d G a ' J ’(r 0 \x 1 ,x a )G t 3 , '{r 0 \x 2 ,x b ) (5.1) 

x E^ iyCTp (r 0 1 x a , x b , x c , x d )G 2 ap 1 s{ r 0 |x c , x 3 , X 4 ) 

Eq. ( 0 ) is a closed integral equation for G 2 , but the 
operator E 2 is given in terms of an infinite series which 
begins with the diagrams shown in FigjKja. 

Equation 0 ) can be considered as a linear inhomo¬ 
geneous equation for G 2 - If we expand around the inho¬ 
mogeneous term G x G we recover the initial expansion 
that was represented diagrammatically above. However, 
we can now seek nonperturbative solutions which ar e th e 
solutions of the homogeneous part of the equation 0 - 
We will demonstrate in subsection A that this nonper¬ 
turbative solution has a power law form in which the 
anomalous exponent A appears, and that this solution is 


di 3 di4G 2 (ro|r,ti = 0,r,i 2 = 0, R, t 3 , R, t 4 )) 

= [r(r,R)G o (ro|r,R,0 + )] 2 g 2 (r,R) ■ (5.4) 

The function g 2 (r, R) is a dimensionless function by con¬ 
struction. In general g 2 (r, R) is a function of r, R and 
the angle between r and R. For R r the angle be¬ 
comes irrelevant. In the regime when r and R are in the 
inertial interval, the function g 2 can depend on the ratio 
r/R only. The notation will be 

y = r/R , g 2 (r, R) = g 2 (y) ■ (5.5) 

In Appendix A we show that this function satisfies the 
approximate integral equation 

52 ( 2 /) = c f —K(y, y a )g 2 {ya) ( 5 . 6 ) 

J y a 
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with some dimensionless constant C, and kernel K(y, y a ) 
which can be written as 


K(y,Va ) 


g(,y/ya)g(,ya ) 

9 ( 9 ) 


(5.7) 


Here 9 ( 9 ) and < 72 (j/) are dimensionless functions which 
are defined by 


g(r/r') = r' 3 G(r Q \r,r',0 + ), (5.8) 

g 2 (r/r') = g 2 {r,r') . 


Eq.( |5.6|) was derived on the basis of the choice of y -C 1 
and considering j <g i/ a < 1 . Of course, the equation 
for the nonlinear Green’s function included additional 
regimes that we did not study explicitly. To gain in¬ 
sight to the r elev ance of these regimes we are going to 
interpret Eq.(5.6) in a more general setting, allowing y 
and y a to go between zero and infinity: 


. 92 ( 9 )=^/ — K(y,y a )g 2 (y a ) • (5.9) 

J 0 ya 


In this subsection we examine the solutions of the model 
Eq.(U). 

We note that the known asymptotic properties of 
G(r 0 |r, r', 0 + ) imply that 


9(2/) = 



for y » 1, 
for y <C 1 . 


(5.10) 


with a and b being dimensionless coefficients. These 
properties imply also the asymptotic properties of the 
kernel K(y,y a ). In the y 1 y a plane we display the evalu¬ 
ation of the kernel in Fig. [2(J. 

In order to understand the type of solutions that are 
supported by the integral equation (|i]) we turn now to 
some simplified models. 


1. Model A 


The simplest model that w e can think of is the one in 
which the constant a in Eq.( 5.10| ) is zero. This leaves 
us with g(y) in the regime y < 1 where our model was 
derived. The asymptotic form of K(y,y a ) is therefore 
KiViVa) =0 for j/o, > 1 and for y a < y. In other words, 
K(y,y a ) = b 2 in region 2 of Fig.^o], and zero in all the 
other regions . W e call this model A. In this case the in¬ 
tegration in ( |5.9|) is the trajectory denoted by the dashed 
line Ai in Fig.| 20 |, which is limited to lie within region 2 . 
The equation for 92 ( 9 ) reads 


92 ( 9 ) = A 0 / — 92 ( 9 a) • (5.11) 

J y y a 


where Ao = Ab 2 . The solution of this equation is 

92 ( 9 ) = G/y A A = A 0 (model A) . (5.12) 


The coefficient C should be determined from the bound¬ 
ary condition 92 ( 1 ). In dimensional form this result reads 

J dt 3 dUG(r 0 \r, h = 0, r,t 2 = 0, R, t 3 , R, f 4 )) 

= C [r(r, R)G°( r 0 |r, R, 0+)] 2 ° . (5.13) 

This is the same result that we obtained in section 3.B on 
the basis of the resummation of the logarithmic contribu¬ 
tions in the simple ladder diagrams. This is not surpris¬ 
ing since model A represents exactly in resummed form 
the nature of the approximation in section 3.B in which 
only the dominant contribution was taken from each dia¬ 
gram. The subdominant contributions are not expected 
to ruin the anomalous scaling behavior, but since the 
anomalous exponent is sensitive to the numerical value 
of the coefficients, we expect the subdominant terms to 
effect the numerical value of the exponent. To study this 
effect we turn to a slightly more complicated model. 


2. Model B 


The next model, which we refer to as model B, is ob¬ 
tained when we use the asymptotic properties of K ( y , y a ) 
as shown in Fig.^O] up to the boundaries of the regions, 
choosing the coefficients a = b. Start with y < 1, and 
follo w th e trajectory Bi of Fig|2(] in the integration in 
Eq. (5.9). This trajectory crosses three regions in Fig]2C|, 
namely 1,2 and 3, and accordingly the equation has three 
integrals: 


92 ( 9 ) = Ao 


1 f y 

— / dy a y a g 2 {y a ) 
y Jo 


rl dy a 


+ I —g 2 {y a )) 
Jy Va 


dy a 

Va 


92 (y a ) 


(5.14) 
, 9 < 1 


The second integral is the one taken into acco unt in 
Model A. We seek a s olutio n in the form of ( 5.12| ). Sub¬ 
stituting this form in (5.14) we find that the requirement 
of convergence of the integrals puts a limit on the allowed 
values of A. These limits are 


—2 < A < 2 . (5.15) 

The two solutions for A are 

A± = 1 ± \/\ — 2Ao (model B , y < 1) . (5.16) 

The two real branches of solutions as a function of Ao 
are shown in Fig.pT|. For A 0 <C 1 the solution A_ of 
model B coincides with the solution o f mod el A. Indeed, 
for small Ao the second integral in ( |5.1S| ) is dominant 
since it is the only one proportional to 1/A. This con¬ 
tribution is equivalent to considering the ladders in their 
fully stretched configuration, which corresponds to model 
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A. However in this case there is another branch of solu¬ 
tions, A_|_, which for small Ao is dominated by the first 
integral in ( |5.14 ) which is proportional to 1/(2 — A). For 
Ao = 1/2 these two branches coincide. In this model 
there are no real solution for Ao > 1/2. 

Since we have two solutions for A, the general solution 
is a sum 


The requirement of convergence is again —2 < A < 2, 
since it is coming from regions 1 and 3 which are the 
same as in model B. Seeking a solution in the form ( 5.12| ) 
we find that A solves the transcendental equation 


j 2 -A 


2-A 


d? 

A 


d A -rf 


2 —A 


1 - A 


1 

Ao 


(5.21) 




C _ 

y A - 


(5.17) 


with coefficients C + and C_ that are determined by the 
boundary condition and the requirement of continuity 
across the boundaries of Fig.26. 


We see that in the regime y < 1 the solution (|5.17| ) 
has a leading scaling exponent which is A+, whereas A_ 
appears only as a correction to scaling . In our simple 
model the leading scaling exponent A + takes on values 
between 2 and 1 before becoming complex. In general we 
may have a whole spectrum of exponents that appear as 
corrections to scaling. 


One obvious property of the solution of this equation is 
that the symmetry A —> (2 — A) is preserved, meaning 
that the solution is still symmetric around the A = 1 line. 
Secondly, in the regime Ao —> 0 we have again the same 
two solutions as in model B i.e. A = Ao and A = 2 — Ao- 
The solution crosses the symmetry line A = 1 just once, 
exactly when 


A 0 = l/2(d + ln d) . (5.22) 

The implication of this result is that the topology of the 
line of solutions remains similar to the solution of model 
B, as shown in Fig.^. In particular the leading anoma¬ 
lous exponent is again larger than 1. 


3. Model C 


B. Exact Solution of the Anomalous Exponent A 


Finally, we discuss for completeness a model C that 
introduces yet more freedom in the scaling function g(y) 


gjy) 

a 


1 

for 

y > l/d, 

Vdy 

for 

d < y < 1/d, 

y 

for 

y <d . 


(5.18) 


with d < 1 (for d = 1 we regain model B). Model C 
introduces a smoo ther cross over in the regime y ~1. 
Substituting (5.18) in the expression for K(y : y a ) du 
we find for y <C d the following expressions 


In this subsection we derive the important exact result 
A = 2 - C 2 ■ (5.23) 


In order to prove this we will first establish an identity 
which relates the mass operator $ which is shown in 
Fig.2b and the two-eddy irreducible mass operator E 2 
which is expanded in Fig.10a. The identity is 


(yq\xi,X2i X3, £4) 


S$ a p(x 1 ,X 2 ) 

SF% C \X 3 ,X 4 ) 


(5.24) 


(: Va/y ) 

[reg. 

1: 

Va<dy 

l> 

dy a /y 

[reg. 

1-2: 

dy < y a 

<y/d], 

1 

[reg. 

2: 

y/d < y a 

< d ], 

d/y a 

[reg. 

2-3: 

d < y a 

< l/d], 

1 /vl 

[reg. 

3: 

l/d < y) 



(5.19) 


The expression for K(y,y a ) in the regions 1,2,3 corre¬ 
sponds to the values in model B, cf. Fig.^. For d = 1 
the regions 1-2 and 2-3 disappear. Now the integral equa¬ 
tion for <72 (j/) takes the form 


g 2 {y) = A 0 
d 

y Jd 


1 

7 

y/d 


r d v 


„2 

y Jo 


dy a y a g 2 (ya) 


dy a g 2 (y a ) 


y 

l/d 


r dye 

fy/d 1J Oi 


g2(y a ) (5.20) 


a r dy “ 1 a _u f°° dy “ / ^ 

+ d / 7^32 iva )) + / —302(1/0) 

J d Va J 1/d Va 


where the superscript “(pc)” means that the functional 
derivative is taken only with respect of the correlators on 
the principal cross section of the diagrams of *I>. 

The proof of the identity is available from inspection of 
Figures 2 and 10. By taking the functional derivative of 
diagram (1) in Fig.2 we get the first contribution to the 
first diagram on the RHS of Fig.10a (The one with two A 
vertices). The factor of 1/2 in Fig 2 disappears because 
there are two identical contributions to the functional 
derivative. If we take all the diagrams in $ with just 
two correlators in the principal cross section, then the 
functional derivative will yield exactly the first term on 
the RHS of Fig.10a. Next consider diagram (2) with 
three wavy lines at the cross section. This is the first in 
an infinite series of diagrams with three wavy lines at the 
cross section. Finding the functional derivatives of these 
contributions we generate all the terms in E 2 with two 
wavy lines at the cross section. For example, the diagram 
(2) in Fig.2 produces three terms, the second, third and 
fourth in Fig.10 with bare vertices. The generalization of 
the procedure is clear. 
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Next we consider the Wyld equation ( |2.9|) , and evalu¬ 
ate the functional 

derivative 5F a p (ro | Xi , x 2 ) /SD^ s c ^ (a ;3 , X 4 ) where the func¬ 
tional derivative has the same meaning as Eq.(2.6), but 
we restrict the variation only to contributions appearing 
in the principal cross section in the diagrammatic repre¬ 
sentation of F. Rewrite the Wyld equation in schematic 
form 


F = G*[D + <S>]*G (5.25) 


where the star operation * means integration over space 
and summation over the tensor indices. In the same 
schematic representation the functional derivative takes 
the form 


8 F 

SD^ “ 
Next observe that 

(5Z)(p c ) 


GG + G* 

<5$ 

— SF (p c ) 


5D(p°) 

(5.26) 

SF 

(5.27) 

* < 5 D(p c ) 


Using now the identity ( 5.24 ) we conclude tha t 
SF/ 8 D^ pc > solves exactly the same integral equation (5.1) 
as the nonlinear Green’s function, i.e the equation dis¬ 
played in Fig. 10b. This means that 


8 F a p(r 0 \x 1 ,x 2 ) 

SD^ s c \x 3 ,X4) 


= G^\r 0 \x 1 ,x 2 ,X3,X4) ■ 


(5.28) 


One should note that if we were evaluating the full func¬ 
tional derivative of F with respect of D , not restricted 
to the principal cross section, we would have derived an¬ 
other relation, which is 

8 F a p( r 0 |a:i,X 2 ) __ / 5 2 w a (r Q \x 1 )w / 3(r Q \x2) \ ^ 

8 Djs{x 3 , 0 : 4 ) \ Sh 1 (x 3 )Shs(x 4 ) / 


Here the RHS is another type of nonlinear Green’s func¬ 
tion. 

At this point we can use the result ( 5.28| ) in Eq.( |5.2| ) 
to obtain the relation 

r fj . 8F a p(r 0 \x-4,x 2 ) <R) 

T a/ 3 (r 1 ,r 2 ,R) = / dx 3 dx A - 7 —t - D jS ( x 3 ,x 4 ) • 

J 8D\ /S '(x 3 ,X4) 


(5.30) 


Now the physical meaning of the function T a p( iq, r 2 , R) 
becomes clear. It is the the change in the double cor¬ 
relator F af 3 (ri,r 2 ) due to the existence of an additional 
random forcing on the scale R with a correlation D 
The point now is that if the coordinates iq and r 2 are 
much smaller than R, then the double correlator should 
have its usual universal exponent () 2 and the role of the 
additional forcing is only in changing the magnitude of 
F. Accordingly 


lim T a/ 3 (r 1 ,r 2 ,R) oc F a ^(r 1 ,r 2 ) ■ (5.31) 

ri,r 2 <Ci? 


Remember that 


F a p(ri,r 2 ) = S a p( ri) + S a p( r 2 ) - 5 a ^(ri - r 2 ) , (5.32) 
S a p(R) = { 8 u a (r + R, r^u^r + R, r)) , (5.33) 


where <5u was defin ed i n (0). Consequently the deriva¬ 
tive im plied in Eq.( |5.3| ) picks up only the last contribu¬ 
tion in ( 5.33 ) with the final result which is Eq. (|5.23 ). 


6. SUMMARY AND DISCUSSION 

In the paper I in this series we presented a proof of lo¬ 
cality in the perturbative theory of the velocity structure 
functions and the velocity Green’s functions. This proof 
excluded the possibility of a perturbative mechanism for 
anomalous scaling behavior of these quantities. In this 
paper we began to explore the nonperturbative origins 
of anomalous scaling in tur bul ence. By examining the 
nonlinear Green’s function fj we showed that its dia¬ 
grammatic series exhibits ladder diagrams that produce 
logarithmic terms that resum to an anomalous power of 
the dimensionless ratio of two separation distances.The 
results of Section 5 can summarized as follows: 

V^G 2 (r,r,R,R) oc r" A . (6.1) 

The analysis of the resummed perturbation series of the 
nonlinear Green’s function in section 5 indicated that A 
has a critical value which is 2 — <^ 2 for which nonpertur¬ 
bative effects may become very important and may lead 
to a renormalization of all the scaling exponents in the 
theory. In Section 6 we demonstrated that A attains ex¬ 
actly this critical value. This result opens up the critical 
scenario for the renormalization of the scaling exponents 
£ n , which will be explored in detail in paper III of this se¬ 
ries. It will be shown there how the critical scenario may 
result in multiscaling with the outer scale of turbulence 
as the renormalization length. The deep implications of 
this type of criticality are explained in 1 22 | and in paper 
III. 
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7. APPENDIX A: REDUCTION OF THE 
INTEGRAL EQUATION 

Consider the homogeneous part of Eq.(|5.l|), integrate 
it over t 3 and £ 4 , and consider it for t\ = t 2 = 0 , 
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r\ = r 2 = r and r$ = = R. Dividing the result by 

[r(r, i?)G°(r 0 |r, R, 0 + )] 2 we find the following equation: 

Mr ' R) = [r(r,B)G°(r„|r,fi,OUl 1 j dx - dndx ^ 

x G(r 0 |r, r a ,t a )G(r 0 \r, r b , t b )T, 2 (r 0 \x a , x b , x c , x d ) 

x J dHdU(hiMx ct x^R,U) ■ (7-1) 

The property of rigidity of £2 implies that the main con¬ 
tribution to the integral comes from the region r a ~ r b ~ 
r c ~ r d . Because of the logarithmic situation the largest 
contribution will come from the region of r a integration 
in which r <C r a R. This offers an evaluation of the 
integrations over r b ,r c and r d as r a 9 . We rewrite 

92{r,R) = 

X J dt a dt b dt c dt d G(r 0 \r,r a ,t a )G(r 0 \r,r a ,tb) (7.2) 

x T, 2 {ro\ra,t a r a ,t b ,r a ,t c ,r a ,t d ) 
x f dt 3 dt<iG 2 (r 0 \ra,t c ,r a ,td, R,h, R,t 4 ) . 

Note now that the typical time scale of G(r 0 |r, r a , t a ) is 
smaller than r(r), whereas the typical time scale of £2 is 
r(r a ) which is much larger than r(r). Accordingly we can 
take t a = t b = 0 in the expression for £ 2 , and pe rform 
the t a and t b integrals with the help of the relation ( [ 2 . 12 ). 
The result is 

92{r,R) = [r(r,i?)G°(r o |r,i?,0+)] 2 / draral1 
x j dt c dt d [T(r, r a )G(r 0 |r, r a ,0 + )] 2 (7.3) 

x £ 2 (ro|r a ,t Q = 0,r a ,t b = 0 , r a , t c , r a , t d ) 

x J dt 3 <lt/iG 2 (ro\r a , t c ,r a ,t d , R,ts, R, t 4 ) . 


To proceed we note that now £2 is determined by one 
time scale, r(r 0 ). On the other hand G 2 depends on 
three time differences, say £3 — t c , f 4 — t c and t d — t c ■ The 
time arguments t c and t d cannot exceed r(r a ) because of 
the decay of £ 2 . For i? 3> r a this means that f 3 — t c « i 3 , 
t 4 — t c « and t d — t c ~ r(r Q ). In the evaluation o f G 2 
we can take t d — t c to be zero. Using the definition (5.4) 
and evaluating the integral over t c and t d as r(r a ) we 
have 


<72 (r, I?) 


i 


dr a r a * l 


[r(r, i?)G°(r 0 |r, R, 0 +)] 2 d 
x [T(r a )r(r,r a )G(ro\r, r a ,0 + )] 2 (7.4) 

x £2 (r 0 1 r a , t a — 0, r a , t b — 0, r a ,t c — 0,r a ,fj 0) 
x [r(r a , R)G°(r 0 \r a , R,0 + )} 2 g 2 (r a , R) ■ 

Using Eqs. ( [2.131 - 2. 16 ) we find that 

r(r a )r(r, r a )r(r a , i?)/r(r, R) = [r(r a )] 2 . 


We need now to evaluate £ 2 .In the previous section we 
showed that all the diagrams in this series are rigid. Con¬ 
sequently the major contribution to the equation (5.1) 
comes from the regime in which the four coordinates in 
£2 are of the same order. Because of this property of 
rigidity and the scaling relations each diagram in the se¬ 
ries has the same order of magnitude and represents a 
homogeneous function of its arguments with the same 
scaling index that we denote as a 2 . We can evaluate £2 
from any one of the diagrams in its series. For exam¬ 
ple the first diagram with two A vertices gives us 1/r 2 
from these vertices, l/(r 6 (r(r) 2 ) from two delta functions 
S(xi — Xj ), and S 2 (r) from the correlator. In total 


£ 2 (r) 


S 2 (r) 


C 


i-(r )] 4 


(7.5) 


where we made use of the scaling relation ( 2 . 15| ), and C 
is a dimensionless coefficient. The last form will be of use 
in our calculations below. The scaling exponent of £2 (r) 
can be read from (|7.5|): 


(7.6) 


o '2 = 2 C 2 — 10 = —6 — 4 z . 


One can check that this scaling relation, together with 
the fact that the scaling exponent of G(r 0 |a; Q , x c ) is -3 
guarantees that the second term in the RHS of Eq. ( 0 ) 
has the same scaling index as t he first term, which is - 6 . 
Using then the evaluation ( 0 ) we find 

f dv ~ 

g 2 (r,R)~A / —K(r,r a , R)g 2 (r a , R) (7.7) 
J r a 


where K (r, r a , R) is the dimensionless function 

" r aG(r 0 |r, r a , 0+)G°(r o |r a , R, 0 + ) n 2 


K (r, r a , R) = 


G°(ro|r, R 1 0+) 


(7.8) 


When all the scales are in the inertial regime I\ (r, r a , R ) 
is a function of two ratios only, 

K(r,r a ,R) = K(y,y a ), y = r/R, y a =r a /R. (7.9) 

In this regime we can write the final equations (^i. 6 |)-(^ 8 |). 
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FIG. 1. Basic graphical symbols for w, BL velocity differ¬ 
ences (2.4); P, vertex - amplitude of interaction; F double 
correlation function of the BL velocity differences (2.7); G, 
Green’s function (2.6) wo, Fq and Go, corresponding values 
in zero order approximation with respect to interaction (“bare 
values”). 


(a) 



(b) 



FIG. 2. Diagrammatic representation of the renormalized 
series expansion for the mass operators, (a) The mass opera¬ 
tor S of the Dyson equation (2.8) and (b) the mass operator 
<f? of the Wyld equation (2.9). 
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FIG. 3. The diagrammatic representation of the B-L ve¬ 
locity. Panel (a) shows the equation for the renormalized 
velocity, and panels (b)-(e) show the expansion up to fourth 
order in the vertex. 
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Fig 4b 
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FIG. 4. The diagrams for the unaveraged response, up to 
fourth order in the vertex. 








FIG. 5. The diagrammatic expansion of the 4-point nonlin¬ 
ear Green’s function (3.1). In this figure we show the simple 
ladder diagrams in panel (a) and the resummation of the sim¬ 
ple ladder diagrams in panel (b) 



FIG. 6. The diagrammatic expansion of the 4-point non¬ 
linear Green’s function (3.1). In this figure we show dia¬ 
grams that add to those shown in Fig.5, and that decorate 
the vertices. The sum of these diagrams turns the vertices 
into dressed ones. 




FIG. 7. The diagrams that dress the vertices. In this paper 
only vertices of type A and B appear. 
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FIG. 8. The diagrammatic expansion of the 4-point nonlin¬ 
ear Green’s function (3.1) In this figure we show the dressed 
simple ladder diagrams, panel (a) and the resummed equa¬ 
tion (panel (b) which results from summing all the diagrams 
in panel (a) . 



FIG. 9. The diagrammatic expansion of the 4-point non¬ 
linear Green’s function (3.1) . Here we show the fully renor¬ 
malized quantity, which is obtained by summing the result 
of the summation in Fig.8 to all the more complex diagrams 
that appear in the series. The dark rungs of the ladder are 
the two-eddy irreducible mass operator, whose diagrammatic 
expansion is shown in Fig.10. 



FIG. 10. (a) the diagrammatic series for the two-eddy irre¬ 
ducible mass operator, (b) The integral equation for the fully 
renormalized nonlinear Green’s function. 



(c) 



FIG. 11. Diagrams used in the explanation of the property 
of rigidity. 
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FIG. 12. Some local geometries about a given vertex that 
appear in the proof of rigidity. 




FIG. 13. Some local geometries about a given vertex that 
appear in the proof of rigidity. 



FIG. 14. The elements that appear in the local geometries 
in Figs. 12 and 13. 


FIG. 15. The final relevant geometries in stretched dia¬ 
grams and the indication where the main contribution to the 
integral comes from. 
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FIG. 16. Diagrams used in the proof of rigidity of the rungs 
of the ladder. 




FIG. 17. Examples of dangerous “ladders within the lad¬ 
der” diagrams. 
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FIG. 18. Ladders drawn in r,t coordinates, to show the 
natural leaning over that is dictated by causality. 


FIG. 21. The solution 5 as a function of Ao in model B and 
model C. 



FIG. 19. More dangerous “ladders within the ladder” dia¬ 
gram. 



FIG. 20. The kernel of models A (trajectory Ai) and B 
(trajectories Bi and B 2 ). 
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